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Standard methods for higher-order calculations of QCD cross sections in hadron- induced 
collisions are time-consuming. The fastNLO project uses multi-dimensional interpolation 
techniques to convert the convolutions of perturbative coefficients with parton distribution 
functions and the strong coupling into simple products. By integrating the perturbative 
coefficients for a given observable with interpolation kernels, fastNLO can store the results 
of the time-consuming folding integrals in tables, which subsequently are used for very fast 
rederivations of the same observable for arbitrary parton distribution functions, different 
scale choices, or a s (Mz)- Various tables with code for their evaluation are available for 
numerous jet measurements at the LHC, the TeVatron, and HERA. FastNLO is used in 
publications of experimental results by the ATLAS, CMS, CDF, DO, and HI collaborations, 
and in all recent PDF analyses by ABM, CTEQ, HERAPDF, MSTW and NNPDF. This 
article focuses on developments implemented in the new version 2 of fastNLO, enhancing 
and broadening its functionality. 

1 Introduction and the fastNLO concept 

Precision measurements in high energy physics reveal their full power only if they are compared 
to accurate theoretical predictions. For the interpretation of experimental data or the extraction 
or tuning of model parameters reasonably fast theory calculations are needed. For measurements 
defined in published cross sections, repeated computations of (almost) the same cross sections 
have to be performed. Typical examples where calculations have to be repeated more often for 
the same observable are: 

• Comparisons of data vs. theory, where e.g. different sets of parton distribution functions 
(PDFs) are used, like those provided by the global fitting collaborations. 

• Determination of the PDF uncertainty on theory predictions. As an example, the un- 
certainty on the cross sections arising from the NNPDF pQ PDF set require 100 to 1000 
rederivations of the same cross sections where only the PDFs differ. 

• Fitting of PDFs where for each step in the iterative procedure the cross sections have to 
be recalculated for the corresponding temporary PDF. FastNLO is used by various global 
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PDF fitting groups, like MSTW [2], CTEQ [3], NNPDF pQ, ABM g], or HERAPDF [5]. 

• Determination of model parameters, or the determination of the strong coupling constant 
in iterative fits. 

• Determination of the theory uncertainty of cross sections from missing higher orders, 
where conventionally the calculations are repeated for different choices of the renormal- 
ization and/or factorization scale. 

• Studies of the scale dependence of the theory cross sections. For processes involving 
multiple scales it is not clear, which scale setting to choose. 

Some observables in next-to-leading order (NLO) or even higher-order calculations can be 
computed rather fast like e.g. DIS structure functions. Other observables like Drell-Yan and 
jet cross sections, however, are very slow to compute. Especially the latter ones are the current 
focus of fastNLO. FastNLO provides computer code and precalculated tables of perturba- 
tive coefficients for various observables in hadron-induced processes. The calculation of the 
fundamental cross sections or matrix elements is performed by flexible computer code like 
NL0Jet++ [51 El [E1IH], which has been used here. 

For illustrating the ideas in this document, all formulae are shown for jet production cross 
sections in deep-inelastic scattering. All concepts and formulae can be generalized to hadron- 
hadron collisions like at the LHC or the TeVatron [ID] . 



Perturbative QCD predictions for observables in hadron-induced processes depend on the strong 
coupling constant a s and on the PDFs of the hadron(s). Any cross section in lepton-hadron or 
hadron-hadron collisions can be written as the multiplication of the strong coupling constant 
to the power of n, a™ , the perturbative coefficients Ci. n for the partonic subprocess i, and the 
corresponding linear combination of PDFs from the one or two hadrons f t , which is a function 
of the fractional hadron momenta x a , x^ carried by the respective partons, as: 



The calculation of this cross section with reasonable statistical precision is very slow due to 
the necessary Monte Carlo integration of the accessible phase space. The idea of fastNLO is 
to separate the PDFs and the a s factors from the perturbative coefficients Ci ; „ and to convert 
this integration into a sum |ll[|12j . This discretization introduces a set of eigenfunctions Ei(x) 
(with ^2 i Ei(x) = 1) around a defined number of x- values. The PDF in eq. [I] can then be 
replaced by fi ~ J^i fa(%i)Ei(x) and is removed from the integral. When calculating the 
perturbative coefficients the nodes receive fractional contributions of each event within the x- 
range of each eigenfunction. The perturbative coefficients are calculated once with very high 
statistical precision and are stored in a table. The remaining integration over x to compute the 
cross section is turned into a sum over the n perturbative orders, i parton flavors, and all the 
x-nodes. To respect the scale dependence, a similar procedure is employed. The multiplication 
of the PDF, a s , and the perturbative factors are performed when reading the table of all Cj jT1 
values, which is very fast and gives the opportunity to change the PDFs and a s as required. 



1.1 The fastNLO concept 
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When calculating hadron-hadron cross sections, the two x integrations are replaced by two 
sums instead of only one like in the depicted DIS case. 

2 New features in fastNLO version 2 

Here, new features of the fastNLO version 2 are presented in comparison to the previous ver- 
sion 1.4 [TU]. The fastNLO project comprises three main elements, the table format, creator 
code to create and reader code to read and evaluate the tables with an interface to PDFs. 
A new table format features more flexibility and foresees to incorporate multiple additive or 
multiplicative contributions to the cross sections. Threshold corrections, which were already 
available previously for hadron-hadron induced inclusive jet production [13] , higher-order calcu- 
lations, electroweak corrections, or new physics contributions can be implemented in a similar 
way as soon as they are available. Further multiplicative correction factors like non-perturbative 
corrections can be stored together with their uncertainties. Also data can be included within 
the new table format together with their correlated and uncorrelated uncertainties. 
The computation of the tables has been optimized. An automated scan to determine the covered 
x-range is performed first. The now flexible number of x-nodes for each analysis bin respects 
their different x-coverage. The scale dependence is stored as a separate array and also the 
interpolation of the scale nodes is optimized. 

The fastNLO reader code for evaluating the fastNLO tables is now available in Fortran as well 
as in independently developed C++ classes with an agreement of 0(1O~ 10 ) between the two. It 
is distributed in one package [TJ] as open source code, which is installable following the GNU 
autotools procedure with the only dependence on some external functionality to access PDFs 
like in LHAPDF [15]. Previously released tables keep their validity and can be converted into 
the new format. 

3 The generalized concept of flexible-scale tables 

A generalized concept for even more flexible tables was also released in fastNLO version 2. 
These are called flexible-scale tables and are based on two principles. 

The dependence on the renormalization and factorization scales ji r / f can be factorized when 
calculating the perturbative coefficients Ci in , like 

Ci,„(Mr, M/) = 4,n + lo g(Mr)c[,„ + log(/i/)c{ n . (2) 

This way, only scale independent weights c°, c r , and c' are stored in three scale independent 
fastNLO tables. The multiplication of the scale dependent log terms are performed only when 
evaluating the table. Similarly to the a s term in eq. [I] where a s (/i r ) can be regarded as an 
arbitrary function of the scale also the scales fi r and /i/ can be regarded as functions of any 
relevant k observables s&, like fJ- r /f = (J-r/ f(s\, S2)- FastNLO examples employ k = 2 since each 
observable needs a separate interpolation array and increases the required evaluation time. 
This method gives the opportunity to store multiple possible scale definitions, like e.g. the jet 
momentum j»t an d the event virtuality Q 2 . When evaluating the fastNLO table it is now pos- 
sible to choose any function of these two scale settings for the definition of the renormalization 
and factorization scale. Typical examples are e.g. fi 2 = (Q 2 + p 2 ? )/2 and fi 2 — Q 2 , but also def- 
initions like /j, = s%- exp(0.3 • S2) are possible. Further, the scales can be varied independently 
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and by arbitrary scale factors. This gives new opportunities to study the scale dependence of 
cross sections. The concept is also available for pp and pp calculations. For future applica- 
tions, the flexible-scale concept is valid also for higher-order calculations like NNLO without 
any significant loss in speed. 

4 Showcase application 

Numerous fastNLO tables are available on the fastNLO website [H] for various measurements 
by ATLAS, CDF, CMS, D0, HI, STAR, and ZEUS. All calculations were performed using the 
NL0Jet++ program [7] [5] 0|] , for calculating the matrix elements. A data/theory comparison 
of global inclusive jet data in hadron induced processes as a function of the transverse jet 
momentum for various center-of-mass energies is shown in fig. [j] employing the MSTW2008 
PDF sets. Hadron-hadron induced processes further include 2-loop threshold corrections, which 
represent a part of NNLO. An updated version of this plot is available through the arXiv 
article [IB] . 
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Vs = 546 GeV 

° CDF cone 0.1 < 

Vs = 630 GeV 

• D0 cone 0.0 < |ri| < 0.5 

Vs = 1 .8 TeV 

o CDF cone 0.1 < 

• D0 cone 0.0 < 

Vs = 1 .96 TeV 

• CDF k T 0.1 < |y] < 0.7 

□ CDF cone 0.1 < |y| < 0.7 

• D0 cone 0.0 < |y| < 0.4 

Vs = 7 TeV 

□ ATLAS anti-k T R=0.6 |y| < 0.3 
t ATLAS anti-k T R=0.4 jyj < 0.3 

• CMS anti-k T R=0.5 |y| < 0.5 

all pQCD calculations using NLOJET++ with fastNLO: 



• ■a' H1 k T Q^: from 150 to 5000 GeV^ 
odao* ZEUS k T Q 2 : from 125 to 5000 GeV 2 



k T Q 2 :from 150 to 5000 GeV 2 
oDAoe ZEUS k T Q 2 : from 125 to 5000 GeV 2 



pp-bar 



PP 

"(X1) 



a s (M z )=0.118 | MSTW2008 PDFs 
NLO plus non-perturbative corrections 
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j pp, pp: incl. threshold corrections (2-loop) 
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Figure 1: Ratios of data and theory for inclusive jet cross section measured in hadron-hadron 
collisions and in deeply inelastic scattering at different center-of-mass energies. The ratios are 
shown as a function of jet transverse momentum px- The theory results are computed for 
MSTW2008 PDFs. 
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